Spin Chirality Fluctuation and Anomalous Hall Effect in Itinerant Ferromagnets 
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The anomalous Hall effect due to the spin chirality order and fluctuation is studied theoretically in 
a Kondo lattice model without the relativistic spin-orbit interaction. Even without the correlations 
of the localized spins, axy can emerge depending on the lattice structure and the spin anisotropy. We 
reveal the condition for this chirality-fluctuation driven mechanism for axy Our semiquantitative 
estimates for a pyrochlore oxide Nd2Mo207 give a finite a^y ~ 10 Q,~^cm~^ together with a high 



resistivity px 



10" 



cm, in agreement with experiments. 
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The quantum charge transport in the background of 
the (fluctuating) spins is one of the central issues in 
strongly corr elated electronic systems This 
issue is relevant to most of the intriguing systems such as 
high-Tc cuprates 0, manganites showing colossal mag- 
netoresistance (CMR) and magnetic semiconductors. 
Of particular interest is the quantum Berry phase gen- 
erated by the solid angle subtended by the spins in the 
noncoplanar configurations, which is called spin chiral- 
ity . Recently it has been proposed that this spin 
chirality contributes to the anomalous Hall effect (AHE) 
in ferromagnet because the spin chirality acts as an ef- 
fective magnetic field or gauge flux b for the charge car- 
riers 10,l|,|l^. These contrast with the conventional 
theories of the AHE, where the spin-orbit coupling is in- 
dispensable. Conventional theories can be roughly cat- 
egorized into two. One is based on the band picture 
where the interband matrix elements between the Bloch 
wave states with the spin-orbit coupling gives the finite 
Hall conductivity axy jl^ . This theory has recently been 
given modern interpretations in terms of Berry phase, 
parity anomaly and Chern numbers [ p^ . The other is 
based on the skew scattering of itinerant electrons by 
fluctuating localized moments 

Experimentally the AHE in Nd2Mo207 has been dis- 
cussed from the viewpoint of the spin chirality p6| , p^ . 
Among other Mo pyrochlore systems, this compound ex- 
hibits a noncoplanar spin structure and a large axy in 
the metallic ferromagnetic state below the Curie temper- 
ature Tq ^ 90 K. axy at a lowest temperature has been 
explained consistently in terms of a mean-fleld theory by 
taking into account the slight tilting of Mo spins coupled 
to the Nd spins However, the temperature depen- 

dence of axy is much slower than that of the transverse 
ordered component M± of the Nd spins, which can not 
be explained in the mean fleld theory pq ]. This experi- 
mental fact, together with the spin-ice nature of the Nd 
spins, strongly suggests that axy is generated by the spin 



chirality fluctuation of the Nd spins in addition to the 
ordered component. 

In this paper we study axy and also the longitudi- 
nal conductivity axx in the presence of the fluctuating 
spin chirality of the localized spins. The purpose of this 
study is twofold. One is to analyze the experiments in 
Nd2Mo207 semiquantitatively. The other is to reveal the 
condition for the spin chirality fluctuation to contribute 
to axy in terms of the lattice structure and the anisotropy 
of the localized spins. 

To be explicit, we start with a model for Nd2Mo207. 
This compound has a pyrochlore structure with two sub- 
lattices of the tetrahedrons of Nd and Mo displaced by 
half a lattice constant. The Nd atoms yield J — 5/2 
localized spins of 4/ electrons, while the Mo atoms two 
itinerant d electrons of the t2g orbitals per site. The lo- 
calized spins are well described by a spin-ice model with 
uniaxial spin anisotropy specified by either parallel or 
antiparallel direction from vertex to center of mass of 
each tetrahedron We treat these localized spins as 
the spin- 1/2 Ising spins. The nearest-neighbor pairs of 
d electrons and localized spins are coupled by an anti- 
ferromagnetic exchange interaction Jfd |l9|] . Then, our 
Hamiltonian reads 

(Xi,X-,),Q Xi 

+ Jfd 51 (1) 

Cxi,Q, 4.^, Tlx. and Sx, = 5 4i,a'^aa'Cx,,Q' are 
the annihilation, creation, number and spin operators of 
conduction electrons on a Mo site with a crystallo- 
graphic index i in the Mo tetrahedron shown in Fig. 
respectively, a is a spin index, a is the vector form of 
the Pauli matrices. Uj is the unit vector pointing to the 
center of mass of the tetrahedron for the Nd spin with a 
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crystallographic index j. S^, = • nj with the spin- 

1/2 locahzed spin operator at a Nd site Xj with an 
index j. - ) denotes the summation over the nearest 
neighbors. The chemical potential /i is self-consistently 
determined with (rii). The itinerant ferromagnetism oc- 
curs due to a double-exchange mechanism in the degener- 
ate Mo 4d orbitals. For simplicity we consider the single 
band model and treat the ferromagnetism in terms of 
the effective magnetic field h^ff in the z direction only 
for the conduction electrons Besides, we assume 

a reduced uniform electron transfer integral t between 
the nearest-neighbor d electrons. In the double exchange 
mechanism [Q ignored here, the coherence temperature 
of conduction electrons is given by Tq. We take into ac- 
count this effect by choosing a reduced value t = 100 
K, which is an order of magnitude smaller than the bare 
transfer integral and is similar to Tq ~ 90 K pl[|. 




FIG. 1. Low-temperature spin structure of the Mo itin- 
erant spins and the Nd localized spins in a unit cell for 
Nd2Mo207. 

We write the Dyson equation for the electron Green's 
function in the 8 x 8-matrix form as G^^{uj + k) = 
Gmf(w -|- i6, k) — S{uj + iS, k) with a positive infinitesi- 
mal number 77, where 

Gr.l,. .,, + ill, k) 

= [(tJ + ^ + Cr^^„,/lcff/2 + ivi)5i^i- + 2iWk,i,i'] 5a,a' 
- JfdSi,i' ^ Wq=0,i J {S^')nj ■ da,a' ■ (2) 
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^<l,i,jXj,j' (Q)'^-q,i' ,i' 

qjj''.a,fi' 

X ((Ta,a • rij) Gi^i>\a,a'{i^ -f i77,k - q) {ffa' ,a' ■ ny) (3) 

fq,i,j is the dimensionless form factor that describes the 
nearest-neighbor couplings between the Mo and Nd sites. 
Xij (q) is the wave- vector dependent susceptibility of the 
localized spins defined by Xjj' (q) = {N /T){5S';i^j5S\y) 
with the temperature T, SS^^^ = S^^j — {Sj')SqQ and 
= TT Ex, 'S'x'^ e-*i-^^ . We note that these locahzed 
Ising spins are static. 

First, we neglect the correlations of the localized 
spins. Then, Eq. (H) is exact with {S^^) = and 



Xj-j' (q) = /4T. Defining the spin components of 
the Green's function and the self-energy part as Gi^i' = 

^2a,a' ^a' ,a.^i,i'\a,a' and Si^i' = X/a,a' ^ a' ,a^i,i' \a-,a' i 

spectively, we obtain 



J 



[uj + irj, k) = ^ Vk-k',i,iW-k+k' 



2(t/i,i' {uj + iri, k') • nj)nj ~ Qi^ {uj + if], k') 



(4) 



Under the itinerant ferromagnetism, Q has a nonvanish- 
ing z component. Since all of rij's deviate from the z 
direction to form a noncoplanar spin structure and every 
Mo site is linked to three pairs of Nd sites with a differ- 
ent index j, the first term in the square bracket in Eq. 

yields a tilting of the conduction electron spins like 
the solid arrows in Fig |l|. This tilting gives a noncopla- 
nar spin structure of the conduction electrons and a spin 
chirality. 

These results can be interpreted as follows. One can 
consider three stages of the approximations, (a) The first 
is the simplest one, namely to take the average of the 
localized spins (Sj) first, and to calculate the spin chi- 
rality for these averaged moments jT^ , which is zero in 
the disordered state, (b) The second is to consider the 
average of the spin chirality for the localized moments, 

-.loc 



(C/i-°^ ) = {Si ■ Sj X Sk) which is again zero for the 



I.e., {un-^} 

disordered case, (c) However, what is relevant to the 
transport properties is the spin chirality of the conduc- 
tion electrons, namely that of Mo spins Si. One should 
then take into account the ferromagnetic moment along 
z direction. Therefore only the x and y components are 
determined by the localized Nd spins. In this case one 
can show that {Ch'i°l^'^') = {si ■ Sj x Sk) does not vanish. 

Now we consider the condition for the applicabil- 
ity of this fluctuation induced a^y First we re- 
place the localized Ising spins by the classical Heisen- 
berg spins. Then, we obtain "Eiji (cj + irj, k) = 



N ^q. 



jGi,i' {uJ + if], k — q) in the local mo- 



ment limit. Here, the tilting effects of the conduction 
electron spins is smeared out, leading to a-xy = 0. This is 
also consistent with the experiment in Gd2Mo207, where 
Gd^+ has no orbital angular moment and hence the Ising 
anisotropy is absent, showing no enhancement of a^y to- 
wards low temperatures p^ ]. Next we consider another 
crystal structure, i.e., the cubic lattice. The Ising spins 
are put on each lattice position pointing toward the cen- 
ter of the cubes. Similar analysis gives again axy = in 
this case. It is because the periodicity of the crystal ex- 
cludes the uniform magnetic field produced by the gauge 
flux, and all the elementary plaquettes are equivalent. 

We can extend these arguments to more general single- 
orbital lattice models under the itinerant ferromagnetism 
with the net magnetization along the z direction, to ob- 
tain the necessary condition as: (I-i) When the transla- 
tional symmetry is not broken, the unit cell of the Bra- 
vais lattice must contain no less than three localized spins 
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with different spin anisotropies and three itinerant spins 
where three are located on different positions for both lo- 
cahzed and itinerant spins. (I-ii) For each itinerant spin, 
the directly coupled localized spins should satisfy the fol- 
lowing two conditions; (I-iia) the spin inversion symmetry 
with respect to the xy plane is broken and (I-iib) there 
remains a finite inplane component of the summation of 
{n^ a^'^^j a)''^] a O'^^^ ^^'^ spins j and the internal degrees 
of freedom of each spin a. For example, in the case of the 
classical spin without anisotropy or with the easy plane 
perpendicular to the z axis, a corresponds to the degrees 
of freedom over (half of) the unit sphere or a circle, re- 
spectively. Then, the integral over this variable gives zero 
and ffxy — 0. (II) If the translational symmetry is bro- 
ken and the unit cell is enlarged, then the new unit cell 
should contain the part satisfying (I-i) and (I-ii). The 
above conditions are also applicable to the case where 
the Ising variables are ordered. 
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FIG. 2. Black lines represent numerical results in the local 
moment limit of the localized spins for (s^), (s^), Oxx and 
cjxy after rescaling [^. Jfd = 80 K and {rn) = 0.7 are 
adopted. Grey lines represent experimental results for Oxy 
and axx at an applied field along the z direction 0.5 T by 
Taguchi et al. jl6[. 



this local-moment case, Uxx, which is proportional to 
[t/ JfdYi is also temperature independent. Moreover, the 
magnitude of axx is similar to experimental results. This 
indicates that the present choice of Jfd/t — 0.8 is almost 
consistent with experiments. 

Next, to make the comparison with experiments in 
Nd2Mo207 more quantitatively, we include correlations 
of the localized spins, which enhance Uxy In the 
present case, the localized spins interact with each other 
through the conduction electrons, as in RKKY interac- 
tion. Xj,j' {^) is calculated by solving the following equa- 
tions: 

Xj^,, (q) = cosh^ ^ + 2JfaV,^,, (q), (5) 

i 

with the ordered component of the localized spin 
determined by (5^^) = itanh(£j72T) with = 
''^q=o,ij (sq=o,i) ■ ^j. The effective interaction 
Vj,ji (q) via the conduction electrons contains the func- 
tion Z\(sq_i), which represents the linear response of the 
itinerant spin with a wave vector q and a crystallographic 
index i to the ordered component of the localized spin 
with a wave vector q and an index j' . Then, calculation 
of Vj,ji (q) corresponds to the sum of the ladder diagrams 
shown in Fig. ^. This term enhances short-ranged corre- 
lations of the localized spins. The short-ranged correla- 
tions reflect fluctuations around the uniform two-in, two- 
out ordered state. We approximately calculate V,j'(q) 
by replacing the momentum-dependent irreducible vertex 
with its summation over the momentum. Namely, mode- 
mode couplings that mix collective excitations with dif- 
ferent momenta are neglected. Then, a newly obtained 
Xi-j'(q) is substituted into Eq. (^) to calculate a new 
Green's function until convergence is reached. 



Therefore Nd2Mo207 is a special case where the chi- 
rality fluctuation of the localized spins contributes to the 
AHE. The fluctuation induced Uxy qualitatively explains 
the experimental results in Nd2Mo207 that cjxy increases 
below Tq even when a ordered moment and spin chiral- 
ity of the Nd localized spins are negligibly small [ p^ . 
We have obtained self-consistent solutions of the Green's 
function and calculated the longitudinal and transverse 
conductivities Oxx and Uxy by neglecting the vertex cor- 
rections. After rescaling the numerical results into the 
system with (n^) = 2, and the J = 5/2 localized spins |2^] 
for comparison with experiments, results for the ferro- 
magnetic polarization of the conduction electrons (s^) 
produced by /icff, amplitude of the inplane spin compo- 
nent of the conduction electrons (s''~), Cxx and axy are 
shown in Fig. ^. Here, with increasing (s^) below Tq — 90 
K, Uxy grows in the absence of correlations of the local- 
ized spins, although it is still two digits smaller than the 
experimental results [l6| . Since the quasiparticle damp- 
ing rate 7 oc J^L/i has no prominent T dependence in 



FIG. 3. Diagrammatic representation of V, .,/(q). Solid 
lines, dashed lines and dots represent the electron Green's 
function Gi^i'ic^^' (k, -I- irf), spin correlation function 
rxjj'(q), and the coupling 2JfdVq^ija ■ UjSi ii. The irre- 
ducible vertex corresponds to the dashed line with two dots 
at the ends. 

Numerical results after the rescaling procedure |2^] 
are shown in Fig. IJ. As (s^) increases below Tc, the 
transverse component of the magnetization M±, (s"*") 
and axy increase. At intermediate temperatures around 
50 K, an enhancement from the results without short- 
ranged correlations is more than 10. There is a crossover 
around ^ Jfd that characterizes a shift from domi- 
nant short-ranged components of the localized spin cor- 
relations at higher temperatures to long-ranged compo- 
nents at lower temperatures. At lower temperatures. 
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axy exhibits a larger enhancement compared with exper- 
iments [0 . These results indicate that inclusion of cor- 
relations, not only the long-ranged order but also short- 
ranged correlation of the localized spins enhance Uxy 
Although the present choice of Jfd/t = 0.4 is half the 
previous one in Fig. a^x remains of the same order of 
magnitude as the experimental results p6| . In a simple 
single-band case coupled to classical spins, 7 increases 
logarithmically with increasing characteristic length of 
spin fluctuations ^. In the present case, molecular fields 
produced by the itinerant spins prevent ^ from diverg- 
ing. Then, below T*, ^ turns to decrease. This yields 
a shallow minimum of Uxx around T^. The obtained 
inferred from the rapid increase of is similar to the 
experimental result |ll6| . Since T* is mainly determined 
from Jfd, Jfd ~ 40 K is also consistent with experiments. 




10 20 30 40 50 60 70 80 90 
T[K] 

FIG. 4. Black lines represent numerical results after rescal- 
ing where correlations of the localized spins are taken ac- 
count. Jfd = 40 K and (n^) =2/3 are adopted. Grey 
lines represent experimental results for a^y and a^x at an ap- 
plied field along the z direction 0.5 T by Taguchi et al. 
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Lastly, we mention relations of the present study to 
other recent theoretical ones. We have shown, in con- 
trast to Ref. ||l^, that there is no uniform gauge flux, 

i.e., /u„itceii(^) ~ ^- Hence the microscopic lattice 
structure is crucial. From this viewpoint, the single- 
band model on a cubic lattice does not show finite Gxy 
even though the second-neighbor hopping is introduced 
in contrast to Ref. . Tatara and Kawamura obtained 
axy proportional to the uniform spin chirality of the lo- 
calized spins within an adiabatic approximation in the 
background of the static localized spins Their ap- 

proximation corresponds to the second stage of approx- 
imations (b) mentioned above. Their theory works in 
the dirty case Jr <^ 1 as in spin-glass systems with an 
exchange interaction J and electron relaxation time t, 
while the present theory works even for Jfd/t < 1 as in 
NdzMozOy 111. 

In summary, we have studied a spin chirality mecha- 
nism of the AHE. We have found that in itinerant fer- 
romagnets with localized spins forming superstructures 
with different uniaxial spin anisotropics, the spin chi- 
rality mechanism works even without the correlations 



of the localized spins. Both long-ranged and short- 
ranged correlations of the localized spins enhance the 
Hall conductivity. Full quantitative explanations of the 
AHE in Nd2Mo207 require including the t2G orbitals and 
Coulomb interactions, which is left for future studies. 
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